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We consider future balloon-borne and ground-based suborbital experiments designed 
to search for inflationary gravitational waves, and investigate the impact of residual fore- 
grounds that remain in the estimated cosmic microwave background maps. This is achieved 
by propagating foreground modelling uncertainties from the component separation, under 
the assumption of a spatially uniform foreground frequency scaling, through to the power 
spectrum estimates, and up to measurement of the tensor to scalar ratio in the parameter 
estimation step. We characterize the error covariance due to subtracted foregrounds, and 
find it to be subdominant compared to instrumental noise and sample variance in our sim- 
ulated data analysis. We model the unsubtracted residual foreground contribution using a 
two-parameter power law and show that marginalization over these foreground parameters is 
effective in accounting for a bias due to excess foreground power at low L We conclude that, 
at least in the suborbital experimental setups we have simulated, foreground errors may be 
modeled and propagated up to parameter estimation with only a slight degradation of the 
target sensitivity of these experiments derived neglecting the presence of the foregrounds. 



1 Introduction 



All models of inflation produce tensor perturbations with an amplitude proportional to the 
Hubble expansion rate during the very early phase of the universe. The ratio of the amplitude 
of tensor to scalar perturbations, denoted by r, is the smoking gun for inflation and an 
important discriminant between classes of inflationary models [1, 2]. 

The possibility of measuring the presence of primordial gravitational waves in the near 
future depends on being able to measure and characterize CMB polarization. The decom- 
position of this field into two orthogonal E and B modes can then be linked directly to the 
properties of the primordial cosmological perturbations. The E-mode polarization pattern is 
generated by all cosmological perturbations while the B-mode is only generated by non-scalar 
sources including gravity waves present at decoupling [3-5]. As a result of this, the detection 
of cosmological B-modes on angular scales on which lensing noise is sub-dominant would lend 
powerful observational support for the presence of inflationary gravitational waves. 

The B-mode signal is a small fraction of the total polarized signal making it especially 
prone to contamination by systematic errors and detector noise. Further complications arise 
from polarized foregrounds which dominate over the B-mode signal across much of the sky [6, 
7], making the detection and characterization of cosmological B-modes a very challenging 
task. Any future claim of a B-mode detection, therefore, has to rigorously show that the 
many potential sources of uncertainty encountered in the CMB data analysis methods are 
not biasing the result. 

The upper limit on the tensor-to-scalar ratio, r < 0.28 (95% confidence level) [8], is 
dominated by the gravitational wave contribution to the temperature power spectrum. Sig- 
nificantly improved constraints on r, however, are only possible with measurements of the 
the B-mode spectrum whose dominant contribution comes from the tensor perturbations. 
There have been a number of studies investigating the feasibility of this measurement in the 
presence of foregrounds and of the analysis techniques that may need to be deployed [9-11]. 

In a recent paper, [12] combined a parametric maximum likelihood foreground cleaning 
method from [13] with a power spectrum estimation method from [14] in a simulated data 
analysis pipeline for ground-based and balloon-borne CMB polarization experiments. They 
found that, depending on the frequency coverage of the experiments, foreground residuals 
might play a role and can potentially bias the recovery of B-modes. In cases in which this 
does not happen, or external information is provided in order to control foreground residuals, 
a Fisher matrix derived estimate indicated that a value of r as low as around 0.04 at 95% 
confidence level might be measured. 

In this work we extend the work of [12] in three important directions: 

1. We incorporate the uncertainties of other cosmological parameters and study their 
effects on our final r estimates. 

2. We study the effect of using a full non-diagonal power spectrum covariance matrix on 
the r determination. 

3. We investigate the possibility of accounting for parameter biases from residual fore- 
grounds at the power spectrum level using a simple two-parameter power law model. 

These extensions allow us to study the impact of foreground cleaning uncertainties and 
parameter degeneracies on the final determination of r. In addition, they allow us to compare 
two different techniques of accounting for the unmodeled residual foreground, templates 
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and priors at the map level and a model of the residual power spectrum at the parameter 
estimation level. 

Throughout this work we assume a ACDM model. The fiducial values of the cosmo- 
logical parameters for our simulations are taken to be that of the WMAP5 best-fit values 
[15]: = 0.0227, Qdmh'^ = 0.111, = 0, h = 0.719, r = 0.084, = 0.963, and 

A = 2.41 X 10~^, where Qb, ^dm and Qk respectively are the density fractions of baryons, 
dark matter and curvature to the critical density, a density needed for the universe to be 
spatially flat, h is the Hubble parameter in units of 100 (km/s)/Mpc, A and rig are the 
primordial power spectrum amplitude and spectral index parameters respectively. For the 
tensor-to-scalar ratio we assume a fiducial value of r = 0.05 , which corresponds to a phys- 
ically interesting goal that several near-term suborbital CMB polarization experiments are 
shooting for. Our simulated surveys are approximately optimized for measurement of this 
level of primordial gravity waves. 

This paper is organized as follows: in Section 2 we review the formalism of measuring 
B-mode polarization and the method we use to clean the raw CMB map from foreground 
contamination. In Section 3 we present our foreground parameter determination methodology 
and results from applying it to simulated data in Section 4, concluding in Section 5. 

2 CMB and foreground polarization 

The aim of this section is to briefly review the main assumptions of this work, starting from 
the physics of CMB polarization and its characterization in terms of 'pure pseudo power 
spectra', then giving a description of our foreground simulations and cleaning technique. 

2.1 Characterization of CMB polarization 

With the onset of decoupling, scalar, vector and tensor cosmological perturbations source 
a quadrupole anisotropy in the photon field. Thomson scattering in the presence of this 
anisotropy then generates CMB polarization. The photons mean free path before recombi- 
nation was too short to allow for generation of any substantial fraction of polarized photons, 
while after recombination, there were few free electrons to act as scatterers. Hence the CMB 
polarization anisotropics we see today mainly originate from the brief period before recom- 
bination. After reionization, an additional polarized signal is produced due to the presence 
of free electrons. In summary CMB polarization depends on the physics of recombination, 
reionization and the primordial perturbations, as embodied in the CAMB CMB anisotropy 
code [16]. 

The CMB field can be fully characterized by using the four Stokes parameters: the 
intensity or temperature T, the Q and U linear polarization vectors and the circular polar- 
ization V. The circular polarization is not produced by Thomson scattering, hence only the 
T, Q and U parameters are used to characterize the CMB. 

The measured Q and U Stokes parameters vary under coordinate transformation as 
headless vectors, and are related to each other by a 45 degree rotation of the coordinate sys- 
tem. The quantity P = (Q + iU), however, transforms as a spin-2 quantity under coordinate 
rotation i.e. {Q±iU) = exp (=p2i'0)(QibzC/), where tp is the rotation angle. The temperature 
is a scalar field and can be expanded using the standard spherical harmonics, Yim- The 
analogous basis for a polarization tensor field is the spin- weighted spherical harmonics basis. 
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Case 


/sky 


Frequency [GHz] 


Noise level [/ui^cMB-arcmin] 


Balloon-borne 


1% 


150, 250, 410 


5.25, 14.0, 140.0 


Ground-based 


2.5% 


90, 150, 220 


10.5, 10.5, 31.5 



Table 1. Balloon-borne and ground-based experimental configurations. All channels have a Gaussian 
beam width of FWHM= 8'. 

Instead of ±2aim-, it is physically suggestive and convenient to introduce their linear combi- 
nations with an even and odd parity nature [3, 5] 

«L = -(-h2a«m + -2«im)/2, (2.2) 

«L = i{+2aim - -2aim)/2. (2.3) 

These coefficients are named E and B modes in analogy with electric and magnetic fields: 
the E-mode pattern is curl free and the B-mode pattern is divergence free. In the absence 
of parity violating interactions, the cross correlation between B and E or B and T vanishes 
because B has the opposite parity of T and E. 

2.2 Pseudo-Cf pow^er spectrum estimation 

The pseudo-Ci power spectrum estimation is one of the most widely used CMB power spec- 
trum estimators. Straightforward application of the this technique to cut-sky polarized CMB 
maps leads, however, to E-to-B leakage. The consequence of this leakage is that cosmologi- 
cally important information contained in the CMB B-modes is overwhelmed by the statistical 
uncertainty of the (much larger) E-modes. 

The pure-pseudo power spectrum estimators, adopted in this work, define a special 
conditions on the boundary of the cut-sky region in such a way that only the pure E and 
B-modes are retained [14, 17-19]. We used the Xpure^ code [14], which is a generalization 
of the XSPECT and Xpol codes [20], for estimating the pure pseudo-power spectra, 
and C^^. Xpure uses the S^HAT library^ [21] - an efficient parallel implementation of 
spin- weighted spherical harmonic transforms. 

2.3 Foreground modelling and cleaning 

In this section we summarize the foreground cleaning of our data analysis pipeline. Further 
details may be found in [12], [13]. Our focus is on typical configurations and sensitivity levels 
of future ground-based and balloon-borne observations, which is summarized in Table 1, and 
correspond roughly to the nominal expected performances of POLARBEAR [22] and EBEX 
[23] experiments. Our analysis method should also be of interest to other experiments like 
SPIDER [24] whose frequency bands roughly overlap with those of POLARBEAR. 

According to previous studies of the intensity of thermal dust emission, the area centered 
around RA = 62° and DEC = —45° is favorable in terms of foreground brightness and 
contrast as well as being accessible by suborbital experiments. This area, therefore, has been 
the target for many past CMB suborbital experiments i.e. Boomerang [25], BICEP [26], SPT 
[27], QUAD [28], ACBAR [29] and QUIET [30]. It is also an area of choice for current and 
future ground-based and balloon-borne experiments, including POLARBEAR and EBEX 
respectively. 

^www . ape . univ-par is7 . fx/ ~radek/pureS2HAT . html 
^www . ape . univ-par is7 . f r/~radek/s2hat . html 
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We consider the frequency interval from 90-410 GHz (see Table 1) and so our Galactic 
sky model contains two diffuse polarized foreground components, namely thermal dust and 
synchrotron. The total intensity simulation is based on observations by [31] and modelling of 
the frequency scalings by [32]. The polarization signal was then modelled as a 10% fraction 
of the total intensity with the pattern of the polarization angles on large angular scales 
given by the WMAP dust polarization template [6]. On smaller angular scales we added a 
Gaussian fluctuation power adopting a recipe first presented in [33]. We refer to [12] for a 
full description of the model. 

As shown in [6] polarization foreground power is comparable or higher than the CMB B- 
mode power on all scales which are accessible to the WMAP experiment. Therefore, a robust 
foreground cleaning procedure needs to be applied to the data prior to B-mode measurement. 

Our approach to this problem is to use the parametric component separation algorithm 
proposed by [13] and applied in [12] for CMB suborbital polarization experiments. This 
algorithm solves the component separation problem in two steps. Firstly, the maximum 
likelihood values of the free parameters of the foreground frequency scalings are determined 
using a spectral index 'profile likelihood' that has been analytically derived under the as- 
sumption of constant foreground scaling across the entire map. Secondly, having fixed the 
foreground frequency scalings, the component amplitudes are estimated in a pixel by pixel, 
using a generalized least-squares algorithm. 

Specifically the spectral index likelihood is given by 



where d is the data, A = A {(3) is a component 'mixing' or frequency scaling matrix with a 
total of Ngpec free 'spectral parameters', /3, which are to be estimated, and N is the noise 
covariance matrix of the data. The maximum likelihood values of (3 are then substituted into 



to solve for the amplitudes, s, of the different components and their covariance Ng- This 
algorithm is implemented with a Markov chain Monte Carlo based maximization in the 
MIRAMARE"^ code. This formalism, like other likelihood based component separation methods, 
has an additional advantage of easily parametrizing and propagating instrumental features 
like calibration uncertainties to the final cleaned maps and covariance matrices. 

As we have seen above, the typical parameters to be estimated in the component sepa- 
ration phase are the foreground and CMB amplitudes, and foreground spectral parameters. 
Given the three frequency channel designs of our balloon-borne and ground-based setups, we 
can only deal with one dominant foreground, which corresponds to the thermal dust for both 
experiments' frequency coverage. In the analysis by [12] it was shown that the contamina- 
tion of the unmodeled and hence uncleaned synchrotron component in the final CMB map is 
negligible for the balloon-borne experiment, owing to the choice of frequency bands. On the 
other hand, for the ground-based case, which has a lower frequency band at 90 GHz where 
the synchrotron emission may still be relevant, the synchrotron contamination is high enough 
to bias the dust spectral index estimate thereby compromising the cosmological B-mode sig- 
nal. Either a prior on the dust spectral index or external information on the synchrotron 

^people . sissa. it/'N^leach/miramare 



2 In Cspec = CONST - (A* N-^dY (A* AT-i A) ^ (A* N'^ d) 



(2.4) 




(2.5) 
(2.6) 
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is therefore required for restoring a cleaning performance comparable to the balloon-borne 
case. 

Although we assume a constant spectral index across the observed region, the likelihood 
Eq. (2.4) can accommodate a spatially varying spectral index, for instance by evaluating the 
likelihood on separate low resolution pixels [34] or by modeling spectral index variation as 
a new component [35]. Both of these foreground model extensions come at the cost of 
introducing of extra map-level foreground free parameters and ideally would benefit from the 
luxury of having extra data channels. Conversely, the approach assumed here be thought of 
as a potentially useful "minimally parametric" approach that may be an appropriate for the 
partial sky, low foreground contrast setting we are investigating. As argued by other authors 
[36, 37] any unmodeled spectral index variations might in the end manifest themselves as 
a residual foreground contamination with a powerlaw spectrum - a theme that we will 
return to in §3.4. 

The main aim of this paper is to demonstrate that instead of relying on external priors for 
the dust spectral scalings or templates for the synchrotron component, one can perform a self- 
contained component separation and then account for diffuse foreground bias in the parameter 
estimation phase through a simple two-parameter residual foreground power spectrum model. 

3 Parameter estimation with foreground cleaned CMB templates 

In this section we discuss the parameter estimation methodology that we have developed for 
the purposes of deriving unbiased constraints on the cosmological parameters in the presence 
of foregrounds. Our methodology incorporates the full power spectrum covariance matrix 
and a model for the residual foreground power spectrum. 

3.1 Likelihood 

The likelihood of the data, compressed into power spectrum bands, given a model is assumed 
to have an offset lognormal distribution [38] 

- 21n£ = J][if ^ - zr]'Q,,[z^^ - z#^], (3.1) 
w 

where zi = ln(C'£), = In(C^), and Qw = C^S^JC^' is the local transformation of the inverse 
covariance matrix, S^J, to the lognormal variables z^. Here Ci and are the estimated and 
model power spectra respectively, and XX denotes the EE and BB spectra. 

For the balloon-borne and ground-based experiments considered here, we used only the 
EE and BB power spectra generated using the Xpure pipeline [14]. The minimum and 
maximum multipoles considered are £min = 20, ^max = 1020. To reduce cut-sky effects, we 
averaged the spectra into multipole bins of width A.i = 40. 

To explore the parameter space of the likelihood, we used the COSMOMC sampler [39] 
in its multi-chain mode. The convergence of the chains was verified by demanding that the 
'Gelman- Rubin statistic' R < 0.1, which measures the variance of the chain means divided 
by the mean of the chain variances. 

3.2 Planck prior 

The suborbital experiments we are considering here are specifically designed to measure the 
EE and BB spectra, which are particularly sensitive to the reionization optical depth, r. 
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and the tensor to scalar ratio, r. However, these experiments are not optimized to make 
CMB temperature measurements and lack the polarization information on the largest an- 
gular scales. We therefore imposed a Planck prior to give parameter constraints that are 
representative of the situation expected when these experiments begin taking or analyzing 
data. More specifically, the joint analysis accounts for the fact that the constraints on r are 
expected to be correlated with the other cosmological parameters, particularly Ug. 

Planck [40], which is currently observing from the L2 Lagrangian point, is expected to 
obtain improved constraints compared to WMAP by a factor of a few on the main cosmo- 
logical parameters. Our simplified Planck TT, TE and EE polarization power spectra, are 
simulated with a modified futurCMB code [41] using the HFI 100, 143, 217 GHz channels 
and with specifications given in [40]. For comparison we also investigated a mock WMAP 
likelihood based on the noise levels of the WMAP Q, V, W bands in temperature and V and 
W bands in polarization [42]. For these likelihoods we used the exact expression from [43]. 

We excluded the satellite BB power spectra from our analysis as our main interest in 
this work is to investigate the constraining power of future suborbital B-mode experiments 
and to develop an appropriate analysis technique. It is likely that a measurement of the 
B-mode from reionization at low £ with Planck will require a more specialized treatment due 
to the higher foreground to signal level [10, 44-46]. 

3.3 Covariance matrix 

The form of the pseudo-C^ covariance matrix in the presence of an instrumental white noise 
and a residual foreground, as defined in Eq. (20) of [12], is given by 

= {2i + i)u, + i2iTT) + ' ^'-'^ 

where the first term in this equation represents the sample plus noise variance of the estimated 
pure B-mode spectra, which are assumed to be diagonal in multipole bin-space. In our anal- 
ysis, this term is computed via 500 Monte-Carlo simulations. The second term corresponds 
to the covariance of the cross correlation between the pure CMB-l-noise and the residual 
foreground power. The C^^^^ is estimated using the model given by Eq. (27) of [12], which 
is where the estimated foreground power spectrum and the uncertainty on the estimated 
foreground spectral indices, A/3, are both factored in. The last term in Eq. (3.2) represents 
the covariance due to the residual foreground, and is fully correlated between different bins. 
This term corresponds roughly to the tensor-product term found in equation (A4) of [13], 
which indeed describes the correction to the noise correlation matrix, Eq. (2.6), due to the 
foreground residuals. Note that this covariance formula in Eq. (3.2) corrects and supersedes 
the one in Eq. (31) of [12]. 

In Figure 1 we show the full covariance matrix computed for the 'balloon-borne basic' 
and 'ground-based template' setups, which are two foreground clean cases investigated by [12], 
and discussed further in the next section. The dominant contribution to these covariances 
are from the variance of the CMB-l-noise spectra, given by the first term in Eq. (3.2). It 
is clear that the off-diagonal elements are at least one order of magnitude smaller than the 
diagonal elements, reflecting the fact that the amount of the residual foregrounds for these 
two setups is subdominant. The overall amplitude of the foreground covariance is controlled 
by a factor and the decay in correlations towards high £ reflects the decaying power 
spectrum of the residual foregrounds. 
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Figure 1. The CMB EE (left) and BB (right) power spectrum covariance matrix for the balloon- 
borne (top panels) and ground-based (bottom panels) experiments in binned-multipole space. The 
colour scale is logarithmic in units of ^K^. The different terms that specify the full covariance matrix 
are given in Eq. (3.2). The off-diagonal elements, which are are several orders of magnitude smaller 
than the diagonal parts, are caused by the presence of correlated residual foregrounds. 



3.4 Modelling the residual foreground 

Foreground residuals in the recovered CMB map resemble the underlying foreground emission 
and bias the estimated CMB signal. If we could reliably account all foregrounds in the 
component separation phase, the resulting effect of the residuals on the estimated CMB 
signal could straightforwardly be propagated in the covariance matrix using Eq. (3.2). For the 
type of experiments we consider here, however, this is not the case since we have insufficient 
frequency coverage to deal with the synchrotron component which hence requires external 
information in the form of a strong prior on the foreground spectral indices or an external 
template. 

The use of templates or external priors on the dust spectral parameters has its own 
limitation as we do not have polarization foreground surveys at frequencies, sensitivities 
and angular scales of interest. Furthermore it is commonly through extrapolation of low 
frequency observations that the required template is produced, but this process is prone to 
astrophysical complications. Ideally we want to perform a self contained analysis. 

Modelling of sources of residual foreground power in temperature at high ^ is of course 
quite common in the literature [27, 47, 48]. We propose that a low (. excess of polarized power 
due to residual foregrounds may be modelled in a similar way using a phenomenologically 
motivated two-parameter power law model. This parametrization is inspired by WMAP 
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modelling [6, 7] and the fact that the dominant contribution to the residual foreground is 
due to the unmodeled synchrotron component, it is given by 

Df^ = AfeL'r(^)'"'' (3-3) 

where Di = i{l + l)C£/27r, Iq is the pivot multipole where the residual model is normalized 
at and Cf^ is the B-mode power spectrum estimate in the lowest bin taken here as a 
normalization constant, ^^fg is the fractional value of the residual amount relative to at 
^0, and /Sj-fg is the residual spectral index parameter. We fit these model parameters together 
with other cosmological parameters in the parameter estimation phase. We normalize D^^^ 
at £o = 40, the central value of the first multipole band in our data. 

4 Results and discussion 

In this section we present our findings concerning the effect of marginalizing over residual dif- 
fuse foregrounds at the parameter estimation step, and concerning the effect of incorporating 
a full non-diagonal covariance matrix on recovered parameter values. 

We considered three cases of map-level component separation. The first two, which 
are representative of self-contained component separation attempts, includes the '66aszc' 
(balloon-borne basic) and ^gbasic^ (ground-based basic) cases. The third case, in which the 
ground-based experiment relies on an external low frequency template to remove the syn- 
chrotron component, is referred to as ^gtemplate^ (ground-based case with template cleaning). 
In Figure 2 we compare the input B-mode power spectrum with the level of residual fore- 
grounds left in the estimated B-mode power spectrum from these three setups. 

We show in Figure 3 the values of the estimated parameters for the above three cases, 
bbasic, gbasic and gtemplate, as well as the Planck only case for comparison. As can be 
seen from this figure, the constraint on r is dominated by the suborbital experiments, while 
the constraints on the other parameters are dominated by the Planck prior. This is further 
illustrated in Figure 4 where we compare the r constraint from bbasic joint analysis and from 
the Planck and WMAP priors alone. The right panel of this figure shows that constraining 
power of the Planck prior suppresses the Ug — r degeneracy, which explicitly validates the 
widely-used single-parameter Fisher matrix r forecasts, also assumed by [12]. 

Returning to Figure 3 it is clear that for both the bbasic and gtemplate cases, the input 
values of all parameters are recovered without bias, while for the gbasic case the estimated r 
value is biased positive by more than 3a. This is because for the balloon-borne or the ground- 
based with template cleaning setups the unaccounted foreground residual is negligible, while 
for the attempt at a self-contained ground-based analysis case, gbasic, the foreground residual 
is significant, and hence adds power to the B-mode power spectrum estimates on large angular 
scales. This extra large-scale power then biases the r value estimate, as shown Figure 7. 

4.1 Foreground control at the level of parameter estimation 

The effect of introducing a residual foreground model Eq. 3.3 in our parameter estimation 
phase is clearly shown in Figure 5 where the previously found bias on r is now eliminated 
and the input values of all the cosmological parameters are recovered with the same precision 
as in the balloon-borne and template cases. In the same plot we show also the gbasic and 
bbasic cases for comparison. The constraints on the residual model parameters are, however, 
poor. 
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Multipole / 

Figure 2. B-mode power spectra of the foreground residuals in tlie CMB map. Tlie ^bbasic' 
(starred blue curve) and ^ gtemplate^ (crossed red curve) cases are for the balloon-borne and ground- 
based (with external template) experiments in which successful foreground cleaning is achieved. The 
^gbasic^ (dotted black curve) case represents a ground-based experiment in which a high level of 
residual synchrotron remains, due to insufficient frequency coverage, leading to parameter biases in r 
if left unmodeled. For comparison we show the model B-mode spectra from tensor modes for r = 0.05 
(solid magenta curve), lensing (dash-dash magenta curve) and the sum of the two (sold cyan curve) 

To understand the effect of incorporating the residual foreground model in our parameter 
estimation, we show in Figure 6 the two dimensional contours of the likelihood. The left 
panel shows the anti-correlation between r and the residual foreground amplitude, A^fg. The 
amplitude parameter reflects the ratio of the residual foreground signal to the total at 
the pivot multipole Iq = 40. From this degeneracy plot we can clearly see the reason why 
the estimated r value from the gbasic is biased: The presence of residual foregrounds at low 
i can mimic a positive tensor-to-scalar ratio. The middle and right panels in Figure 6 show 
the degeneracy between r and /3rfg and ^rfg and /3rfg respectively. In both these cases the 
degeneracies are weak, while the residual spectral index parameter, /3i.fg, is poorly constrained. 
It is therefore safe to assume that despite the poor knowledge of /3rfg, our power law residual 
foreground model can disentangle the residual from the CMB signal and lead to a similar 
precision measurement of r from the self-contained analysis of our simulated ground-based 
experiment. 

To show the above points more clearly, we plot in Figure 7 the data together with the 
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Figure 3. Likelihoods for the seven cosmological parameters from the ground-based and balloon- 
borne experiments. As in Figure 2, the bbasic and gtemplate refer to a balloon-borne and a ground- 
based experiments (with external template) where successful foreground cleaning has been applied. 
The gbasic case shows a parameter bias in the value of r due to foreground residuals. The vertical 
lines correspond to the the input values of the parameters. 



model Ce curves. This ensemble of curves corresponds to models drawn from a converged 
MCMC chain. By comparing the 'tensor' BB curve with the residual foreground, 'rfg', curve, 
we see that the Ci region that is most affected by the residual foreground component coincides 
with the region where most of the gravity wave signal is present, at £ < 200. This plot also 
illustrates the importance of measuring the B-mode power spectrum at I > 200 in order to 
be able to estimate and subtract off the B-mode spectrum from lensing. 
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Figure 4. Effect of priors : imposing Planck like or WMAP like priors on the cosmological parameters 
which are not well measured by the suborbital experiments has negligible effect on the constraint of 
r. This is because the r constraint is mainly dominated by suborbital B-mode data. 




Figure 5. Likelihoods for the ground-based cases: the solid black curve is for the gtemplate, the 
dash-dash red curve is for the gbasic and the dot-dash blue curve is for gbasic with the residual 
foreground model given by Eq. (3.3) included in the parameter estimation. Marginalizing over the 
residual foreground model parameters debiases the gbasic r estimate, r = 0.096 ± 0.019, towards the 
input value of 0.05, obtaining r = 0.064 ± 0.025, which can be compared with r — 0.057 ± 0.021 in 
the gtemplate case. 



4.2 Effect of off-diagonal terms in the covariance matrix 

Up to this point we were using a covariance matrix in which small off-diagonal components 
from the modelled residual foregrounds are ignored. In this section we address the question 
of how incorporating such off-diagonal terms in the covariance matrix affect our results. 

The expression that combines all propagated uncertainties starting from component 
separation to power spectrum estimation is given by Eq. (3.2). The full covariance matrix 
for the bbasic and gtemplate cases is shown in Figure 1. We note that off-diagonal terms in 
both cases are several orders of magnitudes smaller than the diagonal elements. 

In Figure 8 we compare the recovered values of our seven base parameters using diagonal 
{bbasic and gtemplate) and full covariance matrix {bbasic full cov and gtemplate full cov). 
For these cases, the effect of ignoring off-diagonal terms is small, with the shift on r being 
slightly bigger in the gtemplate case. This is because, as shown in Figure 1, the residual 
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Figure 6. 68% ( dark green) and 95% (light green) confidence intervals of the likelihood contours 
illustrating degeneracies among the tensor-to-scalar ratio, r, residual foreground amplitude, A^-fg, and 
spectral index, /3rfg, parameters, for the gbasic case with the residual foreground accounted at the 
parameter estimation level. This illustrates how adding foreground power, Arfg, to the fit debiases 
the r estimate towards the input value of 0.05, and how the uncertainty jS^-fg degrades the statistical 
significance of the r detection. 




50 100 300 600 1000 

Multipole / 



Figure 7. Illustration of the model Cg together with the simulated data for the ground-based ex- 
periment with basic foreground setup, gbasic. The solid lines are the different physical components 
that make up the total model Cg. We found the parameters that maximize the likelihood using an 
MCMC search: the model Ci curves shown here are selected randomly from the MCMC chain after 
initial burn-in. This figure demonstrates how a residual foreground model can be used to subtract 
away excess foreground power at low i. 
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Figure 8. Effect of ignoring off-diagonal elements in the power spectrum covariance matrix on the 
constraints of the cosmological parameters. Plots illustrating the full covariance matrices for both 
the balloon-borne and ground-based experiments are shown in Figure 1. Runs with 'full cov' in the 
legend incorporates a full covariance matrix, while otherwise a diagonal covariance matrix is used. 

foreground level is higher for this case. We conclude that at least for these cases, where 
the off-diagonal terms are small as compared to the diagonal ones, ignoring the off-diagonal 
terms has little impact on our results. 



5 Conclusions 

Foreground cleaning and the propagation of uncertainties related to this process is one of 
the outstanding issues in CMB data analysis for ongoing and future experiments. This is 
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particularly important for CMB polarization measurements and the search for B-modes, for 
which foregrounds are expected to be sizeable at all frequencies and across most of the sky. 

Based on the foreground cleaning procedure presented in [13], we developed and com- 
pleted recent results from [12] concerning the propagation of foreground cleaning uncertainties 
in the main steps of CMB data analysis, with a focus on measuring the tensor to scalar ra- 
tio of cosmological perturbations. We considered typical experimental setups for proposed 
ground-based and balloon-borne experiments, for which the foreground contrast might be 
weak enough to justify the assumption of a spatially constant spectral index built into our 
parametric method. 

These experiments operate at frequencies between 90 and 410 GHz, where the main 
diffuse foreground components are polarized synchrotron and thermal dust emission. We 
examined configurations in which a significant source of bias remains after the foreground 
cleaning phase, due to residual synchrotron emission from lower frequencies. In this case, 
we introduced a method for correcting the data for the presence of such a residual, simply 
by parametrizing the residual foreground power at low ^, and then marginalizing over its 
amplitude and powerlaw index. We found that this technique debiases the r measurement, 
although in the case we studied the foreground parameters are poorly constrained, degrading 
the r detection significance. 

We characterize the error covariance due to subtracted foregrounds, and find it to be 
subdominant compared to instrumental noise and sample variance in our simulated data 
analysis. 

We conclude that for the experimental configurations and template-based simulated 
foreground signals investigated here, it is possible to propagate foreground subtraction un- 
certainties and marginalize over foreground residuals, without biasing the final parameter 
determination results. Therefore, the expected performance of these multi-frequency experi- 
ments, usually quoted in the absence of a properly simulated data analysis procedure, might 
only be marginally degraded by the presence of foreground uncertainties and residuals. 
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